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ABSTRACT 


The  propagation  of  electric  currents  and  voltages  along  a  pair  of 
wires  over  a  ground  plane  is  studied.  The  system  is  assumed  lo  be  non- 
uniform;  i.  e. ,  the  self  and  mutual  inductances  and  capacitances  vary 
along  the  wires.  The  existence  and  uniqueness  of  an  electric  wave  having 
prescribed  initial  values  is  shown  to  follow  from  recent  results  of  the  author 
on  symmetric  hyperbolic  systems  of  partial  differential  equations.  A 
construction  for  this  solution  is  given  in  the  case  of  a  coupler  (1.  e. , 
a  pair  of  wires  which  are  non-uniform  and  coupled  over  a  portion  of  their 
length  only) .  This  coupler  problem  is  reduced  to  a  two-point  boundary 
value  problem,  and  the  latter  is  reduced  to  a  pair  of  Initial  value  problems, 
one  of  which  Involves  a  matrix  Riccati  equation.  A  novel  feature  of  the  work 
is  an  "a  priori"  estimate  which  guarantees  that  a  solution  of  the  (non-linear) 
matrix  Riccati  equation  exists  on  the  whole  line. 
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ELECTRIC  WAVE  PROPAGATION 
ON  NON-UNIFORM  COUPLED  TRANSMISSION  UNES 

Calvin  H.  Wilcox 

Introduction.  This  paper  deals  with  the  propagation  of  electric 
currents  and  voltages  along  systems  of  parallel  wires.  In  the  simplest 
case  there  is  a  pair  of  wires,  or  a  single  wire  over  a  ground  plane.  This 
is  called  a  single  wire  line.  The  case  of  n  pairs  of  wires,  or  n  wires 
over  a  ground  plane,  is  called  an  n-wire  line. 

The  single  wire  line  was  analyzed  by  Kirckhoff  in  1857,  and  the 
associated  "Telegraph  Equations"  have  since  become  the  subject  of  a 
large  literature.  This  case  is  not  considered  here. 

With  the  development  of  more  sophisticated  telegraphy  and 
telephony,  analysis  of  the  multiple  wire  lines  was  needed.  The  steady 
state  solutions  for  uniform  n-wire  lines  have  been  analyzed  by  J.  R. 
Carson  and  R.  S.  Hoyt  [2]  ,  L.  A.  Pipes  [13]  ,  S.  O.  Rice  [17]  and 
others.  The  case  of  non-uniform  multiple  wire  lines  (achieved  by 
tapering  the  wires  or  varying  their  spacing)  has  received  much  less 
attention  ( see  [  7] ,  [  19] ) . 


Sponsored  by  the  Mathematics  Research  Center,  U.  S.  Army,  Madison, 
Wisconsin  under  Contract  No.  DA-1 1-022-ORD- 2059. 
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This  paper  treats  the  propagation  of  transient  and  steady  state 
waves  on  non-uniform  multiple  wire  lines.  For  simplicity,  the  work  is 
presented  for  the  2-wire  case.  However,  this  case  is  typical  of  the 
general  n-wlre  case  and  all  the  results  presented  here  extend  immediately 
to  the  general  case. 

The  propagation  problem  for  non-uniform  2-wire  lines  is  formulated 
in  §1 ,  and  it  is  shown  that  the  existence  and  uniqueness  of  a  solution  is 
guaranteed  by  recent  work  of  the  author  [  21]  under  very  general  conditions. 

The  remainder  of  the  paper  is  devoted  to  a  special  propagation 
problem,  the  "coupler"  problem.  The  2-wlre  line  is  assumed  to  be  non- 
uniform  and  coupled  over  a  finite  portion  of  its  length  (the  coupler)  only, 
and  the  propagation  of  waves  through  the  coupler  is  studied.  The  coupler 
problem  is  formulated  in  §2,  and  a  complete  solution  based  on  the  Laplace 
transform  is  given  in  §§3  and  4.  In  §4,  the  coupler  problem  is  reduced 
to  a  two- point  boundary  value  problem  for  a  system  of  four  first  order 
ordinary  differential  equations.  This  problem  is  then  reduced  to  two 
one -point  (or  initial)  value  problems.  The  reduction  employs  a  matrix 
Rlccatl  equation  satisfied  by  the  reflection  coefficient  matrix. 

Matrix  Rlccati  equations  have  been  studied  by  W.  T.  Reid  [  15,  16]  , 
R.  L.  Sternberg  and  H.  Kaufman  [19]  ,  J.  J.  Levin  [9]  ,  R.  M.  Redheffer  [14] 
and  others.  They  have  been  applied  to  multiple  transmission  line  problems 
by  Sternberg  and  Kaufman  [19]  and  I.  Kay  [7].  Analogous  scaler  Rlccati 
equations  for  the  Impedance  and  reflection  coefficient  of  a  non-uniform 


#376 


-3- 


single  wire  line  were  discovered  by  J.  R.  Pierce  [  12]  and  L.  R.  Walker 
and  N.  Wax  [  20] ;  see  also  [18] .  In  all  of  these  papers,  only  the  local 
solvability  of  the  Riccati  equations  is  shown.  A  novel  feature  of  the  present 
work  is  an  "a  priori"  estimate  which  guarantees  the  existence  of  the  desired 
solution  of  the  Riccati  equation  on  every  interval. 

The  special  case  of  a  "directional  coupler"  ,  which  is  of  great 
interest  for  applications  [10]  ,  is  analyzed  in  §5.  Finally,  in  §6  the 
analysis  is  adapted  to  the  steady  state  case. 
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S 1 .  The  Propagation  Problem  for  Non-Uniform  Coupled  Transmission 
lines.  The  electric  currents  and  voltages  on  a  pair  of  coupled  loss-less 
transmission  lines  are  governed  by  equations  of  the  form  [  1 ,  13] 


+  L 


31, 


1  at 


m  8t 


!!i 

8x 


=  0 


f 


"2  at 


8i. 


+  L 


m  8t 


®®2 

8x 


=  0 
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(1-  1)  / 


8e  ae,  ai. 
^llt"  ^m'8r''*’8x 


0  , 


ae  ae.  ai 
^2“  ‘^m‘at”'‘’'a3r  °  • 


Here  t  is  the  time,  x  is  a  coordinate  aldng  the  lines,  and  1^  and  i^ 

(resp.  e^  and  e^)  are  the  electric  currents  (resp.  voltages)  on  the 

first  and  second  lines,  while  and  are  their  self  inductance  per 

unit  length,  C,  and  C,  their  capacitance  per  unit  length,  and  L  and 

their  mutual  inductance  and  capacitance  per  unit  length.  The  lines 

are  assumed  to  be  non-uniform;  1.  e.  L,,  L,  ,  L  ,  C, ,  C-  and  C  are 

1  2  m  1  2  m 


functions  of  x. 
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The  quadratic  form 

"  ■  i  'Vl‘2  +  V2  *  °A  *  «=mV2  +  °2*1  ' 

is  assumed  to  be  positive  definite  and  is  interpreted  as  an  energy  density 
(energy  per  unit  length)  on  the  lines  [13]  .  The  quadratic  form 

is  interpreted  as  a  "Poyntlng  vector"  describing  the  flow  of  energy  (energy 
per  unit  time  in  the  x-dlrectlon)  on  the  lines.  The  equations  (1.  1)  imply 
the  conservation  law 


in 

8t 


+ 


8L 

8x 


=  0  . 


The  most  general  systems  of  linear  partial  differential  equations 
that  possess  quadratic  energy  densities  and  Poyntlng  vectors  related  by  a 
conservation  law  were  described  by  K.  O.  Friedrichs  [  6]  who  called  them 
"symmetric  hyperbolic"  systems.  The  system  (1.  1)  assumes  the  symmetric 
hyperbolic  form  when  written  in  the  matrix  notation 
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The  propagation  problem  for  coupled  transmission  lines  is  to  find  a 
solution  of  (1.1)  having  prescribed  initial  currents  and  voltages, 


i^(x,  0)  =a^(x)  ,  i^Cx,  0)  ej{x,  0)  =bj(x),  e^Cx,  0)  =b2(x) 

for  -00  <  X  <  oo  .  Local  existence  and  uniqueness  theorems  for  the  initial 
value  problem  for  symmetric  hyperbolic  systems  have  been  given  by  K.  O. 
Friedrichs  [6],  P.  D.  Lax  [8],  G.  F.  D.  Duff  [5]  and  others.  A  global 
existence  theorem  for  solutions  with  finite  energy  was  given  by  R.  S. 
Phillips  [11].  Recently,  the  author  [21]  has  proved  a  global  existence 
and  uniqueness  theorem  for  solutions  with  locally  finite  energy;  i.  e. 
satisfying 


f  q(x,  t)dx<«o  for  all  finite  x. ,  x  . 
’‘l 


(1.2) 
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The  initial  values  are  required  to  satisfy  (1.  2)  with  t  s  0 ,  and  the 
coefficients  are  required  to  satisfy  the  conditions 


(1.3)  Lj{x),  L2(x),  L^(x),  C^(x),  C^Cx)  and  C^(x)  are  bounded, 
Lebesgue  measurable  functions  of  x  on  -<»<  x  <  «> ,  and 


(1.4)  there  exists  a  positive  constant  6  such  that 


for  all  ip  i^,  e^,  and  all  x  on  -oo  <  x  <  «  . 
Note  that 


m,  .2 


m  ,  ,2 


m  ,2 


=  h'*!  *  rf  ‘2''  ♦  ¥*2  *  r‘  ‘l'-  *  °1<*1  *  cf  »2'‘  *  '=2<'2  * 

*  i  X  z 


j2  i2  22 


'2  “1 


2  '"1 


Hence  if 


— 

L  L  - 
12 


<k" 


^1^2  - 


(1.5) 
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than 


Thus  (1.4)  holds  provided  L^,  L^,  Cy  are  positive  and  bounded 
away  from  zero  and  (1. 5)  holds  with  a  constant  k  <  1  . 

If  the  Laplace  transform  with  respect  to  t  is  applied  to  (1.1) 
a  system  of  four  linear  first  order  ordinary  differential  aquations  is 
obtained.  The  solution  with  locally  finite  energy  guaranteed  by  [  21] 
can  be  constructed  in  terms  of  the  solutions  of  this  system.  This 
method  is  used  §§3  and  4,  below,  to  solve  the  propagation  problem 
in  a  special  case,  the  coupler  problem. 
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§2.  Formulation  of  the  Coupler  Problem.  In  the  remainder  of  this 

report  the  transmission  lines  are  assumed  to  be  coupled  on  a  finite  segment 

0  <  X  <  i  only,  and  to  be  uniform  and  uncoupled  outside  this  segment,  in 

the  sense  that  L  =C  =0  and  L,,  L, ,  C,  and  C,  are  constants  for 
mm  1  6  1  c 

x<0  and  x>i  (the  same  constants  on  both  intervals).  Under  these 
circumstances  the  segment  0  <  x  <  i  will  be  called  a  "coupler"  between 
the  two  lines.  It  can  be  shown  that  energy  can  pass  from  one  line  to  the 
other  only  in  the  coupler. 

The  voltages  and  currents  on  the  uncoupled  portions  of  the  lln^s 
satisfy  the  simple  transmission  line  equations 


(2,1)  ( 


=  0 


f 


8e  8i 
c  — ^  4*  ^  s  0 

^4  As.  •  ^ 


) 


8x 


where  ,  i^  and  e^  refer  either  to  the  first  line  (j  =  1)  or  the  second 

line  (J  =2).  This  system  is  equivalent  to  the  system 

"  °J8x^^^jS  ■  ®j)  =  0  » 
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where 


and 


are  the  chairacterlstlc  "wave  velocity"  and  "impedance"  of  the  lines, 
respectively.  It  follows  that 


(2.  2) 


I  ij  =f(cy  -  X)  +  o{Cjt  +  X)  , 

< 

e^  =Zj(f(Cjt-x)  -o(Cjt  +  x))  , 

L 


where 


(2.3)  f(c^t-x)  =  Z^i^  +  ej/2Z^  and  g{c^t  +  x)  =  Z.i^  -e^/2Zj 

are  uniquely  determined  by  the  currents  and  voltages  in  the  lines. 

If  hj(T)  a  h^CT)  e  0  f<  r  t<0  then 


(2.4)  (i^,  i^,  e^,  e^)  =(hj(Cjt-x),  h2(C2t-x),  Zjhj(Cj^t  -  x) ,  Z2h2(C2t-x)) 

solves  the  system  (1.1)  for  -oo  <  x  <  « ,  t  <  0  .  This  solution  may  be 


interpreted  as  a  wave  propagating  along  the  portions  x  <  0  of  the  lines 
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toward  the  coupler  and  striking  It  at  a  time  >  0 .  The  propagation  of 
the  wave  through  the  coupler  is  described  by  the  solution  of  the  Initial 
value  problem  for  (1.1)  with 


(2.  5)  {lj(x,  0),  0),  ej(x,  0),  e^{x,  0))  =  (hj(-x),  h2(-x),  Zjhj(-x),  Z^h^C-x)) 


This  problem  is  called  "the  coupler  problem"  in  what  follows. 

The  current  and  voltage  on  line  1  will  have  the  form  ( 2.  2)  for 
x<0,  -oo<t<<»  where,  by  (2.3)  and  (2.4),  f(c^t-x)  =h^(Cjt-x) 

and  g(c^t +x)  =  0  for  x<0,  t<0.  Thus  f(T)  s  hj(T)  for  all  t 

and  g(T)  s  0  for  t<0.  Similar  reasoning  applies  to  line  2  for 

X  <  0  and  to  lines  1  and  2  for  x  >  i  .  It  follows  that  the  solution 
to  the  coupler  problem  has  the  form 

(ij>  ^2’  ®r  ®2^ 

(h^(Cjt-x)  +  rj(Cjt  +  X) ,  h2(c2t-x)  +  r2(C2t  +  x) ,  Z^{hj(Cjt-x)  -rj(Cjt  +  x)}  ,  Z2{h2(c2t-x)  -r2{C2t  +  x)}  ) 


for  X  <  0  ,  -00  <  t  <  00  ,  where  r,(T)  s  r  (t)  a  o  for  t  <  0  ,  and 

(^1>  ^2’ ®1’ ®2^ 

for  X  >  i  ,  -00  <  t  <  00  ,  where  tj(T)  s  t2(T)  =0  for  t  <  -i  .  Evidently, 
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the  functions  r^(c^t  -f  x)  and  ■¥  x)  may  be  Interpreted  as  reflected 

waves  produced  by  the  coupler  on  lines  1  and  2,  respectively,  while 
t^(c^t  -  x)  and  t^Cc^t  -  x)  may  be  interpreted  as  waves  transmitted  through 
the  coupler. 

The  solution  to  the  coupler  problem  may  be  derived  from  the  solutions 
for  two  special  cases  by  a  superposition  principle  (Duhamel's  Principle;  cf. 
[4,  p.  202]).  To  describe  it  let  H(t)  denote  Heaviside's  function,  i.  e. 


H(t) 


1, 

0, 


T  >  0  , 
T  <  0  , 


and  let  (i|,  e|,  e^)  be  the  solution  of  the  coupler  problem  corresponding 

to  an  incident  wave  ( 2. 4)  with 


(2.6) 


hj(T)»H(T),  h^CT)  a  0 


2  2  2  2 

and  (1^,  i^,  e^,  e^)  the  solution  corresponding  to  an  incident  wave  (2.4) 
with 


(2.7) 


hj(T)«0,  h2(T)aH(T). 


Then  It  is  easy  to  verify  that 
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( i^(x,  t) ,  i^Cx,  t) ,  t) ,  e^Cx,  t)) 

00 

(2.8)  »  ^  /  (i[(x,t-T),  i^Cx^t-T),  e}(x,t-T),  e2(x,t-T))hj(CjT)dT 

00 

+  ^  /  (lj(x,t-T),  l2(x,t-T),  ej(x,t-T),  e2(x,t-T))h2(C2T)dT 


defines  the  solution  of  the  coupler  problem  for  an  arbitrary  Incident  wave 
(2.4).  Hence,  the  remainder  of  this  report  treats  the  cases  (2.6)  and 
(2.7)  only. 
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S3.  Formulation  of  the  Laplace  Transformed  Coupler  Problem.  In 
what  follows,  the  Laplace  transforms  of  ti  e  solutions  (l!,  iLe|,  e^)  and 

1  M  1  6 

6^,62)  are  denoted  by  (ij,  I*,  eJ,  E^)  and  (I^,I^,Ej,E2)j  l.e., 

00 

I,(x»8)  =  /  e"®*lj{x,t)dt  ,  etc. 

0 

The  solutions  to  the  coupler  problem  have  the  form 

,,1  ,111. 

=  (H(Cjt.x)  +  j{Cjt  +  x),  r2j(C2t  +  x) ,  Zj{H(c^  t-x) -r^ ^(c^t  +  x)}  ,  -Z2r2j(C2t  +  x)), 


2  2  2  Z 

t^l»^2’®l’®2^ 

=  »  H(C2t-x)  +  r22(C2t  +  x) ,  +  x),Z^  {H(C2t-x)  +  x)} ) 

£orx<0,  where  *  ^22^"^^  “  °  for  x  <  0  .  Hence, 


application  of  the  formula 
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(3.1) 


i  ^  00  T 

f  e‘®*f(ct*  x)dt  =  l/c  e  °  /  e°  f(T)dT 
0 


gives 


sx 


sx 


sx 


sx 


a  -  (1 +e  ^  Rjj(s),  e  ^R^jCs),  Z^{l  -  e  Rj^{s)},-Z2e 


(3.2) 


(if.I^.Ef.E^) 

sx  sx  sx  sx 

=  1  (e  ^  R^2<s)»  1  +  ®  ^R22(®)» ’^1®  *^22^®^^^ 


for  all  X  <  0 ,  where 


ST 

V»  J  « ■ 


Similarly,  on  x  >  i 
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“Z'ZZ'^Z' 


where  tj^(T)  a  a  »  0  for  T<-i  , 

whence  by  ( 3. 1) 


,l  I  , 

=-  (e 


s(x-i) 
c. 


'I 


s(x-i) 

c. 


8(X-i) 

c. 


Tjj(s),  e 


T2i(s),Zje 


Tji(s),  Z^e 


T2i(s)) 


(3.3) 


,,2  .2  _2  _2,  1  , 

=  g  (e 


c. 


S(x-i) 


S(x-i) 


'1 


.1^ 


e 


T22(s),  Zje 


Tj2(s).  Z^e 


T22(s)) 


for  all  x  >  i  ,  where 


s  e 


si 


ST 

00  "  c 
/  e 
-i 


j 


t,|^(T)dT 


I  <  i,k  <  2  . 


On  0<x<i  the  behavior  of  (l|,  E^)  *  J  =  2  ,  is  governed 


by  the  Laplace  transform  of  the  system  (1.1).  Making  use  of  the  well-known 
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rule 


/  ^  dt  =-£(0)  +  8/e’®^f(t)dt 


and  the  initial  values  (ij,!^, ejje^)  »  (i^, l^jOj, e^)  “(OjO 
tsO,  0<x<i,  gives 


5r*»'¥l*V2'''> 


(3.  4) 


j  j 

dr-"®<Vi"W  =« 


‘^^i  J  J 

,0 


dl 


V 


^  +s(C^e|4C,E^,)  =0 


,0,0)  for 


for  0  £  X  <  i  and  J  =  1,  2  . 

Any  locally  Integrable  solution  of  ( 3.  4)  will  satisfy 
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e|(x)  =e}(Xq)  -  s/  {Lj(e)i}(e)  +  6)1^2* ^)}de 

*0 

where  the  intearal  is  finite,  by  (1.3),  and  similar  equations  hold  for 
E^2*  l|  and  *  Hence,  any  such  solution  is  (absolutely)  continuous. 
In  particular,  (3.2)  and  (3.3)  imply 

(i}(0),I2(0),e}(0),E2(0)) 

=  7(1  +  Rll(s),  R^its),  Zj{1-Rjj(8)}  ,  -Z2R2i(8)) 

(3.5) 

(lf(0),l2(0),Ef(0),E2(0)) 

=  i(Rj2(s),l  +  R22(s),  -ZjRj2(s).Z2{1-R22(s)})  , 


and 

(l}(i),l2(i),E}(f),E2<^)) 

fa'®'-  ■ 

(3.6) 

(lf(i),l2(l),E^(i),E2(i)) 

=  'j(Tj2(s),  122(8),  ZjTj2(s),  Z2T22{s)) 


#376 


-19- 


It  Is  shown  below  that  the  differential  equations  (3.4)  and  the  boundary 
conditions  (3.5)  and  (3.  6)  determine  uniquely  the  Laplace  transforms  of 
the  currents  and  voltages  within  the  coupler,  together  with  the  reflection 
coefficients  Rjj^(s)  and  transmission  coefficients  Tjj^(s).  Moreover, 
it  is  shown  that  the  Rj^(s)  can  be  determined  directly,  without  determining 
Tjj^(s),  e|(x,  s),  etc.  Notice  that  once  Rj]^(8)  is  known  the  currents  and 
voltages  in  the  coupler  can  be  obtained  from  the  solution  of  the  initial  value 
problem  (3.4),  (3.5),  and  the  values  at  x  =  i  of  this  solution  determine 


T^,^(s),  by  (3.6). 
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§4.  Solution  of  the  Laplace  T  ransformed  Coupler  Problem.  The 


matrix  notation 


1 

'  *1 

u 

II 

,  E  = 

1 

'V 

is  used  in  what  follows.  Equations  (  3.  2)  may  be  then  written  in  the  matrix 


form 


/ 


\ 


Ri2(s) 


R22(s) 


/ 


sx 


e 


/ 


22 


(s) j 


or  simply 


#376 


-21- 


R(s)  is  called  the  reflection  coefficient  matrix.  Similarly,  equations  (3.3) 
may  *^0  written 


-22- 
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where 


/I„C)  I,2C)\ 


T(8) 


VT2i(s) 


^22^ 


may  be  called  the  transmission  coefficient  matrix. 

The  differential  equations  ( 3.  4)  also  take  a  simple  form  in  matrix 
notation,  namely 


(4.3) 


where 


sU  =  0 


< 

^  +  sCE  =  0 
dx 


for  0  <  X  <  1 


and 

C  = 

Cm\ 

Vl 

^  m 

hi 

\c 

'  m 

^2/ 

The  continuity  conditions  (  3.  5)  and  (  3.  6)  take  the  form 
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(4.4) 


and 


I  sl(0)  =  U  +  R(s) 
sE(0)  =  ZIU  -  R(s)]  , 


(4.5) 


sl{i)  =T(s)  , 

sE(i)  =2T(s)  . 


Equations  (4.  3),  (4.4)  and  (4.5)  do  not  constitute  a  proper 
boundary  value  problem  for  I(x)  and  E(x)  ,  as  they  stand,  because 
R(s)  and  T(s)  are  not  known.  However,  (4.4)  and  (4.5)  imply 


(4.6) 


ZsI(O)  *  sE(0)  =  2Z 

< 

Zsl(f)  -  sE(i)  =  0 

k 


It  will  be  shown  that  (4.3)  and  (4.  6)  define  unique  functions  I(x), 
E(x). 

4> 

If  M  is  an  m  X  n  matrix,  let  M  denote  the  complex  conjugate 
of  the  transposed  matrix.  Note  that  L  =  L  and  C  =  C ,  because  L 
and  C  are  real  and  symmetric.  It  follows  that  if  E  and  I  satisfy  (4.3) 


then 
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(4.7)  J 


+  sI*L  =  0  , 

dx 


^  +  sE*C  =  0  . 
dx 


#376 


This  Implies  the  following 

Lemma  4. 1  If  E  and  I  are  2X1  matrix  solutions  of  (4.3)  on 
a  <  X  <  b  then 

(4.8)  E*I  +  I*EI^  =-2(Res)  /  (  I*LI  +  E*CE)dx  . 

O 

a 


Proof.  By  (4.3)  and  (4.7) 


f 


E  +  I 


dx 


=  -  sI*LI  -  sE*CE  -  sE*CE  -  Sir'Ll 
=  -  2(Res)(I*LI  +  E*CE)  . 


If  z  =  X  +  iy  Is  a  2X1  matrix  with  real  and  imaginary  parts 


X  and  y ,  respectively,  then  a  simple  computation  gives 
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and 


£ 

=  ,  V)V  i  “ 


2*cj=  i;  >0 . 

j  )  K-i 


Combining  this  with  the  preceding  lemma  gives 

Corollary  4. 1  If  E  and  I  are  2X1  matrix  solutions  of  (4.3) 
and  a  <  X  <  b  and  Re  s  ^  0  then 

(4.9)  (E*I  +  I*E)^^  <  (E*I  +  I*E)j^^^  . 

This  result  Implies 

Theorem  4. 1  ( Uniqueness  Theorem) .  The  2 -point  boundary  value 
problem  (4.  3),  (4.6)  has  at  most  one  solution,  provided  Res  >  0. 

Proof.  Suppose  that  (E',!')  and  (E",  I")  are  two  solutions.  Then 
E  =  E'  -  E"  and  1  =  1'-  I"  will  solve  (4.  3)  and 


(4.10)  E(0)=-ZI(0),  E(i)=ZI(i). 


Let 
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= 


/e{\ 

V  V 


h[\ 


V  V 


,  j  =  1,  2 


Then  E^.  are  2X1  matrix  solutions  of  (4.3)  and  satisfy 


(4.11)  E^(0)  = -ZI^(O)  ,  E^(i)=ZI^i),  j=l,2, 

by  (4.10).  Applying  (4.9)  to  E^  and  on  0<x<j  and  substituting 
( 4. 1 1 )  gives 

2(I^*(l)ZI^(i)  +  I^*(0)ZlV))  1  0  ,  j  =  1,  2  . 

But  Z  =  ((6jj^Z^))  with  Zj  >  0  and  Z^>0,  from  which  it  follows  that  1^(0)  =0 

and  therefore  also  E^(0)  =0,  by  (4.  11)  .  This,  with  (4.3),  Implies  that 
E(x)  =  I(x)  =  0  on  0  <  X  <  J  ,  by  the  uniqueness  theorem  for  the  initial  value 
problem  for  (4.3).  Thus  E'  =E",  I'  =  I"  on  0<x<f  ,  which  proves  the 
uniqueness  theorem. 

The  form  of  the  boundary  conditions  (4.  6)  suggests  a  change  to  the 
new  dependent  variables 

E  =  |(I  +  Z'^E),  1  =  |(I-Z“^E) 


which  satisfy 
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E  =  jZ(E  -  f),  I  =  j(E  +  f)  . 

A  simple  calculation  gives 

Lemma  4.  Z  E  and  1  solve  the  2 -point  boundary  value  problem 
(4.  3)t  (4.6)  if  and  only  if  t  and  f  solve  the  2 -point  boundary  value 
problem 

^  +  s(AE  -  BI)  =  0 

or  0  <  X  <  1  > 

^  +  s(BE  -  AI)  =  0 
dx 

where 

A  =  4  (CZ  +  Z’^L)  ,  B  =  4  (CZ  -  Z'^L)  , 

and 


(4.13)  E(0)  =  U  ,  f(i)  =  0  . 

Moreover, 

f(0)=R(s).  E(i)=T(s). 


The  2- point  boundary  value  problem  (4.12),  (4.13)  can  be  solved 
by  solving  two  1-point  (initial  value)  problems.  This  reduction  is  based  on 
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Lemma  4. 3  If  E  and  1  solve  ( 4. 12)  and  exists  on  an 
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interval  a  <  x  <  b  then 

(4.14)  R  «  R(x,  s)  «  f(x)fi(x)”^ 
solves  the  matrix  Riccati  equation 

(4.15)  ~  -t-sdaR- AR-  RA-I-B)  -0 

on  a  <  X  <  b  . 

Proof.  Differentiating  the  matrix  equation  AA  ^  s  U  gives  the  rule 

-1,  -1  ,  -1  —1 
dA  /dx  a -A  dA/dxA  .  Differentiating  the  equation  RsiE  and  using  this 

rule  and  (4.12)  gives  (4.15). 

A  converse  to  this  result  is  provided  by 

Lemma  4. 4  If  R  solves  (4.15)  and  E  solves 

A 

HP  ^ 

(4.16)  g  4  S(A-  BR)E  a  0 

then  E  and  I  s  RE  solves  ( 4. 12) . 

Proof.  Putting  RE  «  f  in  (4.15)  gives  the  first  of  equations  (4. 12). 
Computing  df /dx  =  dR/dx  E  4  RdE/dx,  from  (4.15)  and  (4.16)  gives  the 


second. 
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The  function  R(Xt  8)  defined  by  (4.14)  satisfies  the  boundary 
condition  R(i,  s)  =0,  by  (4.13).  Thus,  the  solution  to  the  2-polnt 
boundary  value  problem  (4.12),  (4.13)  can  be  constructed  by,  first, 
Integrating  (4.15)  from  x  =  i  to  x  =  0  starting  with  R(i,s)  =0  and, 
second,  Integrating  (4.16)  from  x«0  to  xsj,  starting  with  E(0)  sU. 
This  construction  will  fall  only  If  the  function  R(x,  s)  defined  by  (4. 15) 
and  R(  i,  s)  =0  falls  to  exist  on  the  whole  Interval  0  <  x  <  i  .  It  Is  well 
known  that  solutions  of  non-llnear  equations  may  have  this  behavior.  However, 
In  the  present  case  It  Is  ruled  out  by  the  following 

Theorem  4.  2  If  Re  s  >  0  then  the  solution  R(x)  of  (4.15) 
satisfying  R(  I )  s  0  satisfies  the  Inequalities 


(4.17) 


'  ZjIRjj(x)|Sz2IR2i(x)|^  <  7.J  , 

zJr^2<’‘>>'  +  2JR22(x)I^  <  Z^  . 


on  any  Interval  5  5  x  <  i  on  which  It  exists. 

Proof.  Define  i(x)  on  &5x<i  by  (4. 16)  and  i(6)  =  U  | 
and  put  f(x)  =  R(x)E(x).  Then  E,  f  solve  (4.12),  by  Lemma  4.4,  and 
hence  E  =  Z(E  -  I)  and  I  =  E  +  I  solve  (4-.  3)  on  |  <  x  <  J  and  have 


the  values 
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(4.18)  I(6)»U  +  R(e),  E(e)  «Z{U  -  R(e)) 

and 


I(i)  =E(i),  E(i)  =ZE(i)  whence  E(i)  >ZI(i)  . 
In  component  form,  equations  (4.18)  read 


m  = 


/l  +  Rll(S) 


^  R2i(4) 


R^^U) 

1  +  R22(e) 


and 


/z^(l-Rjj(|))  -ZjRj2(e)  \ 


E(|)  = 


V 


Now,  the  two  column  vectors  E^  and  I^  (j  =1,  2)  from  E  and  I  satisfy 
(4.  3)  on  4  <  X  <  1  .  Hence,  by  corollary  4. 1  , 


(4.19)  (E^*I^ +  <  (E^*I^ +  I^*E^)^^g  ,  j=l,2. 

Taking  j  =  1  and  substituting  the  above  boundary  values  gives 
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2I^*(i)ZlNi)  <Zj(l 


Rj  1  (I ))(!  +  Rj  1  (6 ))  -Z2(R2i<e ))  +  (1  +  Rj  1  (6 ))  Zjd  - Rj  J  (I )) 


.  2[2j-ZjIRjj({)I^  -  22IR2j(6)I^]  . 


This  gives  the  first  of  inequalities  ( 4. 17 )  ,  since  Z  is  positive  definite.  The 
second  is  derived  similarly  from  (4. 19)  with  j  =  2  . 

The  reduction  of  the  2-potnt  boundary  value  problem  (4.12),  (4.13)  to 
two  Initial  value  problems  is  guaranteed  by 


Theorem  4.  3  The  initial  value  problem 


(4.  20) 


'  dR 

~+s(RBR-AR-RA  +  B)=0  ,  0<x<i  , 

R(i)  =0 


has  a  unique  solution  R  =  R(x,  s)  for  every  s  with  Res  >  0  .  The  initial 
value  problem 


(4.  21) 


^  +  s(A-  BR)E  =  0,  0<x<l  , 

E(0)  =U 
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has  a  unique  solutlont  and  E  and  f  «  RE  provide  the  (unique)  solution  to 
the  2-polnt  boundary  problem  ( 4. 12) ,  ( 4. 13) . 

Proof.  The  urrlqueness  of  R(x,  s)  follows  from  a  standard  theorem 
on  the  Initial  value  problem  for  differential  equations  [  e.  g. ,  3,  Ch.  2] . 

The  existence  of  R(x,  s)  on  0  <  x  <  i  follows  from  Theorem  4.  2  and  the 
standard  existence  and  continuation  theorems  [  3,  Ch.  2] . 

Corollary  4.  2  The  reflection  coefficient  matrix  R(  s)  *  R(  0,  s) 
can  be  obtained  directly  from  the  solution  of  the  initial  value  problem  (4.  20). 

A  A 

E(  x)  and  I(  x)  can  then  be  obtained  as  the  solution  of  the  initial  value 
problem  for  (4.12)  with  E(0)  » U,  1(0)  -R{a)  . 

Corollary  4. 3  ZR(s)  »  R*(  s)Z  . 

Proof.  ZR(s)  sZR(0,  s)  and  XsZR(x,  s)  satisfies  the  equation 
^  +  s(XHX  -  KX  -  XK*  +  ZHZ)  =  0 

H  =BZ’^  *  • 

K  «ZAZ'^  »  •|(ZC  +  LZ"S  . 

Evidently,  H*  s  h  and  therefore  X*(x,  7)  satisfies  the  same  equation  as 
X(x,s).  Since  X*(i,  s)  ■X(i,8)  sO,  It  follows  that  X*(x,  7)  aX(x,8) 


where 


and 
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or  0  <  X  <  i  ,  by  the  uniqueness  theorem  for  the  Initial  value  problem. 

Corollary  4.4  * 

Proof.  Corollary  4.3  Implies  that  *^2^21^^^’  *'^0'^®ver, 

It  Is  easy  to  verify  that  R^^i  analytic  for  Res  >  0  and  is  real  for  real  s  , 
which  Implies  that  ^21^^^  °  ^21^  ‘ 

The  reflection  coefficient  matrix  can  be  obtained  by  solving  (4.  20) 
numerically.  For  this  purpose,  the  matrix  Riccatl  equation  should  be  written 
as  a  System  of  four  ordinary  differential  equations.  On  carrying  out  the  matrix 
operations  this  system  is  found  to  have  the  form 


*^*Hl  2 

djT  ®^®ll^l''‘®12*4l’^2l‘*’®21^2*Hl'*'®22^2*^2r^l^r^2‘^2r^l^r^21^2‘'‘®ll^  *  ° 


/ 

(4.  22) 


^^^2  2 

dx  ®^®ll’Sl^2  ■‘■®12‘Hi‘^2.  ■‘■®21^2 ■*■^22^2*^22  ‘  ^1^2  "  ^2*^22  '  '^2^1  ‘  ^22^2  ''■^12^ 


‘^'^21  2 

dx  ®  ^®11*^21*Hi  ®12*^21  ®21*^22*Hi  ®22*^22*^21  ’  ^21*Hl  "  ^22^21  ‘  ^1*^21  ’  ^21*^22  ''■^21^  “  ° 


dx  ^*^11*^21^2  ^12*^21*^22  ® 21*^22*H2  ''■^22*^22  "  ^21*^2  ’  ^22*^22  '  ^2*^21  ‘  ^22*^22 ''' ® 22*  “  ° 

V 
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5  5 .  Directional  Couplsf .  The  "directional"  couplers  are  of 
particular  interest  for  applications  [  10] .  They  are  characterized  by  the 
property  that  the  reflected  wave  le  zero  on  line  1  ( line  2)  and  the  transmitted 
wave  Is  zero  on  line  2  ( line  1)  whenever  the  Incident  wave  is  zero  on  line  1 
(line  2)  .  An  equivalent  property  is  ■  0  and  t^^iT)  ■  t2j(T)  ■  0 

for  all  T .  The  corresponding  property  of  the  reflection  and  transmission 
coefficient  matrices  is  *  ^22^^^  *  ^12^^^  *^21^*^  *  ^  for  all  s. 

It  is  known  that  sufficient  conditions  for  a  coupler  to  be  directional  are 

[10] 

(5.1)  L =  -c /vc;^  on  0  <  X  <  i  , 

mit  mi6 

and 

( 5.  2)  “  ^2  on  0  <  X  <  1  , 


where  and  Z^  are  the  constant  Impedances  on  the  uncoupled  portions  of 
the  lines.  The  sufficiency  of  these  conditions  is  verified  below. 

The  matrices  A  and  B  can  be  written 


*-2 


V 


1  “I 


ly"  +  Z.C  Z.C. 

Z^  1  m  2  2 
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In  this  case  the  matrix  Rlccatl  system  (4.22)  becomes 
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dT  *  >'-zf  ■  f,  =  » 


“"iz 


^  jf  Vzz-z;  4  -  +  '^2'^2  -  - 


=  0  , 


(5.3)/ 


dR,,  L  ,  L  _  L 

-*^'-^''1i-z7‘'22»11-<'^1+'^2>"21-zJ>  •“  • 


dx 


dR  L  L  _ 

♦=<-17  V22-z7"22"i2-^'^2'>22'  =“ 


ax 


These  equations  imply 

Theorem  5. 1  If  (.5.1)  and  (5.2)  hold  then  Rjj(s)  =R22(s)  =  0, 
R^^{s)  =  {Z^/Z^)R^^{s)  ,  and  R^^t  s)  =  •^12<  is  the  value  at  x  =  0  of  the 

solution  = 


dR, 


*r  -  "12  + 1 '^1  + '^2"'u  +  - 


0  <  X  <  f  , 


R^2(f.s)  =0  . 

\ 


(5.4)  { 
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Proof^  If  Is  defined  by  (5.4),  s)  is  defined 

to  be  (Zj/Z^) 

zero  then  the  resulting  matrix  R(x,  s)  is  easily  verified  to  satisfy  (5.  3) 
and  R(i,s)  =0.  But  these  conditions  define  R(x,  s)  uniquely,  whence 
R(x,  s)  must  have  the  form  Indicated. 

Corollary  5. 1  If  (5.1)  and  (5.2)  hold  then  the  coupler  is 
directional  and  the  reflection  coefficient  Rj^(s)  is  characterized  by 
(5.4). 


Rj2(x,  s),  and  Rj^(x,  s)  and 


R22(x,s) 


are  defined  to  be 
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§6.  Formulation  and  Solution  of  the  Steady  State  Coupler  Problem. 
The  steady  state  solutions  of  the  transmission  line  equations  (1. 1)  have 
the  form 


yx,t)  =  Re{e^‘^  Ij(x)}  ,  e^ix.t)  =  Re{  e^*^  Ej(x)  }  ,  J  =  1,  2  , 
where  E^,  E^  satisfy  the  equations 


dE 

ar  W  -0 


f 


dE. 


^T—  +  ML  I,  +  L,I,)  =  0 
dx  ml  2  2 


f 


^+MC^E^  +  C/,)  =0 


^  +  MC^E^  +  C,E,)  =0 
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In  terms  of  the  column  vectors 


/e  \ 

,  I  = 

*2 

\  ) 

v  / 

the  equations  take  the  more  concise  form 


(6.1)  < 


dE 

^  +  1«L1  =  0  , 

dx 


d  I 

2i+l«CE=0  . 


In  discussing  steady  state  problems  it  is  convenient  to  deal  directly  with  the 
complex- valued  vectors  E(x)  ,  I(x)  rather  than  the  real -valued  vectors 
e(x,  t)  =Re{e  E(x)}  ,  i(x,  t)  =  Re  {e  I(x)}  . 

On  uniform  uncoupled  portions  of  the  transmission  lines,  the  steady 
state  solutions  have  the  form 


(6.2)  { 


^  -ikjX  ikjX  ^  -ikjX  ikjX 

=  a^e  +  b^e  ,  E^  =  Z^(  a^e  -  b^e  )  , 


+‘’2®  ’  ^z'Z^fa^e  -b^e  )  , 
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where  a^,  a^,  and  are  constants,  and  =  «/c^  ,  k^  =  w/c^  . 

In  the  case  of  a  coupler  occupying  the  segment  0  <  x  <  i  ,  the  constants 
a,,  a,,  b,,  b,  may  have  different  values  on  the  two  uncoupled  portions 
X  <  0  and  x  >  i  . 

The  steady  state  coupler  problem  Is  to  find  a  colutlon  of  ( 6. 1) 

corresponding  to  a  prescribed  steady  state  wave  incident  on  the  coupler 

from  the  left.  This  means  that  for  x  <  0 ,  a^  and  a^  have  prescribed 

value  (and  and  b^  are  to  be  determined),  while  for  x  >  i  ,  b^  =  b^  =  0 

( and  a^  and  a^  are  to  be  determined) .  The  solution  for  arbitrary  prescribed 

values  (a^  a^)  on  x  <  0  can  be  derived  from  the  special  cases  ( a^,  a^)  =  (1,  0) 

and  (a^,  a^)  =  (0,1)  by  superposition.  Hence,  only  these  cases  are  considered 
.^1  “^2  ”^2 

here.  If  E  ,  I  and  E  ,  I  denote  the  solutions  for  these  two  cases 
then  for  x  <  0 , 
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while  for  x  >  i  , 


^1  -lkj(x-i) 
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In  matrix  notation  these  take  the  form 
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for  X  <  0  and 


forx>i.  In  particular,  I  =  I(x,w)  and  E  =  E(x,w)  satisfy  the 
differential  equations  ( 6. 1)  for  0  5  *  ^  ^  boundary  conditions 


(6.3) 


and 


(6.4) 


I(i)  =Q(«)  , 

E'(i)  =ZQ(«)  , 


where 
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P(«)  - 

and  Q(w)  ^ 

are  (unknown)  reflection  and  transmlsaion  coefficient  matrices. 

Problem  (6.1),  (6.3),  (6.4)  is  formally  Identical  with  problem 
(4.3),  (4.4),  (4.5)  of  §4  if  s  =  lu  and 

(6.5)  I(x,  w)  slul(x,  iu),  E(x,  w)  -iwE(x,  iu)  . 

Moreover,  the  second  problem  was  shown  to  have  a  unique  solution  when 
Re s  >  0  ( Theorems  4. 1  and  4.?),  This  proves 

Theorem  6. 1  The  steady  state  coupler  problem  has  a  unique 
solution  which  is  given  by  ( 6.  5) .  Moreover, 

P(«)=R(iw),  Q(w)  =T(iu>) 

and  therefore  P(w)  is  the  solution  of  the  initial  value 

'  dP 

^  +  1«(  PBP  -  AP-PA+B)»0,  0<x<i  , 

P(i,«)  =  0  . 


-44- 


#376 


Note  that  equality  holds  In  (4.9)  when  RbsbO.  Hence,  (4.9) 
implies  ( see  the  proof  of  Theorem  4.  2) 

Corollari/  6.1  The  reflection  and  transmission  coefficient  matrices 
for  the  steady  state  coupler  problem  satisfy  the  identities 

Zj(|Pjj(«)l^+  IQjj(v>)1^)  +Z2(|P2i(«)l^+  IQ2i(«)I^  “^1  ' 

Zj(  I P^  2(«)  I  ^  +  1 0,  2<")  1  ^)  +  Z2(  I  P22<“)  I  ^  +  I  ^22^"^  •^>=22  • 
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